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Abstract 

A  set  of  stereological  procedures  has  been  developed  for  the  rigorous 
quantification  of  microstructural  features  resolvable  using  scanning  electron  microscopy 
in  a+p  processed  a/p  titanium  alloys.  This  paper  identifies  the  four  microstructural 
features  that  most  likely  influence  the  mechanical  properties  in  a+p  processed  titanium 
alloy,  including:  the  size  of  the  equiaxed  alpha,  the  volume  fraction  of  the  equiaxed 
alpha,  the  volume  fraction  of  total  alpha,  and  the  thickness  of  the  Widmanstatten  alpha 
laths.  The  details  regarding  the  quantification  methodologies  are  provided,  as  are  the 
origins  of  the  associated  uncertainties. 

Keywords:  microstructural  quantification,  stereology,  titanium 

1.  Introduction 

The  balance  of  properties,  particularly  the  high  strength,  good  ductility,  and 
enhanced  fracture  toughness,  of  a+P-processed  Ti  alloys  is  particularly  attractive  in  many 
applications  in  the  aerospace  industry.  It  is  claimed  that  these  properties  are  dependent 
upon  the  microstructural  features  [1-5].  However,  there  does  not  exist  a  model  which 
accurately  predicts  the  mechanical  properties  of  a  microstructure  given  a  set  of 
microstructural  variables.  Indeed,  the  very  best  tools  available  establish  whether  a 
microstructural  feature  has  a  positive  (+),  negative  (-),  or  neutral  (0)  influence  on  the 
mechanical  properties  of  Ti  alloys  [5].  The  absence  of  any  physically  based  model  is 
largely  predicted  by  the  complicated  interdependencies  that  exist  among  the 
microstructural  features.  A  solution  to  the  complex  problem  of  isolating  and  quantifying 
the  magnitude  of  the  functional  dependencies  of  property  on  microstructure  or  processing 
involves  the  use  of  rules-based  approaches.  One  such  approach  is  neural  networks  based 
on  Bayesian  statistics  (i.e.,  Bayesian  Neural  Networks)  based  on  the  work  of  D.  MacKay 
[6-8].  This  approach  is  particularly  powerful  in  that  the  neural  network  is  optimized 
without  over-fitting  for  a  particular  problem.  The  fidelity  of  a  solution  (e.g.,  predicted 
property)  provided  by  any  rules-based  approach  is  dictated  by  the  quality  and  amount  of 
the  data  used  to  train  it.  Hence,  for  the  development  of  neural  networks  based  on 
databases  relating  microstructures  and  mechanical  properties,  it  is  essential  to  develop 
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rigorous  procedures  for  providing  an  accurate,  quantified  description  of  the 
microstructural  features  present  in  a  digital  micrograph. 

Previously,  the  authors  have  presented  such  rigorous  stereological  tools  in  a 
newly  developed  computer  imaging  analysis  toolkit  for  P-processed  a/p  Ti-alloy  type 
microstructures  [9-11].  They  have  demonstrated  the  possibility  of  developing  automated 
or  semi-automated  tools  to  determine  average  values  of  the  complex  microstructural 
features  present  in  these  alloys,  such  as  the  thickness  of  the  Widmanstatten  a-laths,  the 
colony  scale  factor,  the  prior  P  grain  factor,  and  the  volume  fraction  of  Widmanstatten  a. 
The  same  authors  later  revisited  the  importance  of  detennining  the  volume  fraction  of 
colony  in  a  microstructure  which  contains  both  basketweave  and  colonies,  and  as  a  result, 
developed  a  microstructure-property  neural  network  based  model  with  a  predictive 
accuracy  of  ±  3%  of  an  experimentally  measured  property  -  a  value  that  has  defined  the 
state-of-the-art  modeling  of  properties  for  Ti-based  alloys  [12].  Indeed,  the  stereological 
procedures  developed  for  p-processed  a/p  Ti-alloy  type  microstructures  has  been 
demonstrated  to  have  overcome  many  of  the  limitations  of  other  manual  procedures,  such 
as  user  subjectivity  [11]. 

The  goal  of  this  paper  is  to  extend  the  original  work  to  include  a  set  of  tools  to 
describe  quantitatively  the  microstructural  features  present  in  a+p  processed  a/p  Ti- 
based  alloys,  such  as  Ti-6-4.  While  some  of  the  features  of  the  a+p-processed 
microstructures  are  similar  to  those  of  the  p-processed  microstructures  (e.g., 
Widmanstatten  a-laths  described  by  an  average  thickness),  others  are  quite  different,  as 
seen  when  comparing  the  p-processed  microstructures  (Fig.  1(a))  with  the  a+p-processed 
microstructures  (Fig.  1(b)).  The  most  distinctive  microstructural  feature  uniquely  present 
in  the  a+p  processed  material  is  equiaxed  a  (often  called  globular  a  or  primary  a).  Such 
equiaxed  a  forms  during  the  thermomechanical  processing  and  subsequent  aging  in  the 
a+P  phase  field,  and  is  reported  to  play  a  crucial  role  in  dictating  the  mechanical 
properties  for  a  given  microstructure.  Indeed,  to  a  first  approximation,  the  presence  of 
such  a  has  been  argued  to  lead  to  enhanced  fatigue  properties  due  to  its  ability  to  blunt  or 
arrest  cracks  [4]. 

Based  upon  both  the  previous  work  describing  p-processed  Ti  alloys  (for  those 
microstructural  features  common  to  both  p-processed  and  a+p-processed  Ti  alloys)  and 
the  current  qualitative  understanding  of  deformation  mechanisms  in  a+P-processed  a/p 
Ti  alloys  such  as  Ti-6A1-4V,  the  microstructural  features  that  have  been  identified  for 
quantification  are: 

1 .  mean  equiaxed-a  size 

2.  volume  fraction  equiaxed-ot 

3.  volume  fraction  total  ot  and, 

4.  thickness  of  Widmanstatten  a-laths  in  transformed  p. 

The  current  paper  describes  the  development  of  the  stereological  procedures  for 
these  four  microstructural  features.  Although  the  estimation  of  the  thickness  of  the 
Widmanstatten  a-laths  has  been  described  in  previous  work,  it  is  included  here  for 
completeness,  with  an  additional  discussion  regarding  the  limitations  to  the  original 
derivation. 
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2.  Experimental  procedures 

A  total  of  34  samples  (2  samples  from  17  unique  thennal  histories)  from  one 
industrially  processed  heat,  and  hence  of  a  uniform  composition,  were 
thennomechanically  processed  in  the  a+p  phase  field.  The  composition,  nominally  Ti- 
6A1-4V,  was  determined  to  be  Ti-6.5Al-4.0V-0.15Fe-0.190  using  inert  gas  atomization. 
The  billets  were  forged  using  industrial  practices  at  968°C  to  two  strain  levels  (s  =  40% 

&  80%),  sectioned  into  samples,  and  subsequently  subjected  to  combinations  of  different 
solution  temperatures  (>  968°C),  cooling  rates  (0.37°C/sec  to  3.7°C/sec),  aging 
temperatures  (593°C  to  704°C)  and  aging  times  (2  to  8  hours).  Table  1  provides 
additional  details  regarding  the  intermediate  thermomechanical  processing  conditions. 

The  samples  were  subjected  to  room  temperature  tensile  tests.  Following  the 
tensile  tests,  an  undeformed  section  of  the  grip  from  each  sample  was  sectioned  for 
metallographic  preparation.  Following  traditional  metallographic  techniques,  the  samples 
were  characterized  using  a  FEI/Philips  Sirion  scanning  electron  microscope  (SEM) 
operating  in  backscattered  mode  at  15  kV  with  a  resolution  of  approximately  2.0  mn. 

Each  image  had  a  resolution  of  3872  x  2904  pixels  and  a  depth  of  8  bit.  The  resolution 
afforded  using  the  SEM,  especially  when  compared  with  other  available  techniques  (e.g., 
optical  micrographs),  is  of  paramount  importance  to  provide  the  highest  fidelity 
quantified  microstructural  data  to  the  neural  networks.  The  specimens  were  imaged  at 
four  random  locations  avoiding  overlap  and  edge  effects.  These  four  micrographs  were 
quantified  using  the  developed  stereological  techniques  (described  below)  to  improve 
sampling  statistics  and  to  provide  an  overall  representation  of  the  specimens’ 
microstructures.  The  micrographs  were  saved  using  the  uncompressed  tagged  image  file 
format  (*.tiff). 

The  backscattered  electron  micrograph  in  Fig.  2(a)  clearly  shows  the  two  distinct 
phases,  namely  the  a  (hep)  and  [!  (bee)  phases,  which  are  most  easily  identified  by  their 
relative  contrasts  -  a  result  of  their  respective  average  atomic  mass.  The  a  phase, 
stabilized  by  elements  with  a  low  atomic  number  (e.g.,  A1  and  O)  appears  darker  relative 
to  the  p  phase,  stabilized  by  the  heavier  elements  (e.g.,  V  and  Fe).  In  addition  to  the  Z- 
contrast,  there  is  also  channeling  contrast  that  is  present,  and  is  most  easily  observed 
within  the  a  phase  ( see  fig.  2a).  Figs.  2  (b-c)  are  histograms  of  the  contrast  variation  in 
the  alpha  phase  with  and  without  contrast  variations  due  to  channeling  contrast, 
respectively,  while  Fig.  2  (d)  shows  the  complete  histogram  for  the  microstructure  in  Fig. 
2  (a).  Given  the  nearly  continual  contrast  variation  that  is  observed  in  these 
microstructures,  and  the  fact  that  the  stereological  procedures  are  developed  based  on  a 
clear  distinction  between  the  two  phases,  it  is  essential  to  threshold  the  micrographs  to 
precisely  delineate  the  two  phases.  Therefore,  using  a  semi-automated  thresholding 
procedure,  including  a  Gaussian  operation  to  reduce  the  channeling  contrast  present  in 
the  a-phase,  each  micrograph  was  converted  to  a  binary  image  (see  Fig.  2  (e)). 

Additional  preprocessing  steps  are  required  for  each  of  the  stereology  procedures 
which  follow.  Of  paramount  importance  is  the  proper  identification  and  selection  of  the 
equiaxed  a.  Such  a  task  presents  difficulties  as  the  equiaxed  a  that  is  present  in  these 
alloys  can  exhibit  a  significant  range  in  size,  depending  on  the  thermomechanical 
processing,  as  shown  by  some  of  the  extreme  variations  in  size  and  volume  fraction  of 
microstructural  features  that  can  occur  in  Figs.  3(a-c).  These  three  figures  correspond  to 
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samples  forged  within  the  a+P  two  phase  field  at  968°C  to  either  40%  strain  (Figs. 

3(a,c))  or  80%  strain  (Fig.  3(b)).  Following  the  forging,  the  sample  corresponding  to  Fig. 
3(a)  was  solutionized  in  the  two  phase  field  at  985°C,  cooled  at  3.7°C,  and  aged  at  704°C 
for  two  hours.  The  sample  corresponding  to  Fig.  3(b)  was  solutionized  in  the  two  phase 
field  at  977°C,  cooled  at  0.37°C,  and  aged  at  593°C  for  eight  hours.  The  sample 
corresponding  to  Fig.  3(c)  was  solutionized  in  the  two  phase  field  at  968°C,  cooled  at 
0.37°C,  and  aged  at  704°C  for  two  hours.  The  microstructural  variation  that  is  illustrated 
in  Figs.  3(a-c)  presents  a  significant  challenge  to  the  development  of  an  automated 
routine  for  the  identification  and  isolation  of  the  equiaxed  a.  Currently,  the  most  robust, 
automated  routine  will  either  miss  small  equiaxed  a  particles  or  identify  as  equiaxed  a 
regions  that  would  otherwise  be  considered  to  be  grain  boundary  a  or  a  large 
Widmanstatten  a-lath.  An  example  of  this  misidentification  is  shown  in  Fig.  4  (b)  for  the 
representative  micrograph  shown  in  Fig.  4(a).  Therefore,  additional  manual  selection  of 
equiaxed  a  is  required  to  define  all  equiaxed  a  particles  (Fig.  4(c))  and  then  separate  the 
individual  particles  (Fig.  4(d)).  Work  is  currently  underway  to  develop  a  fully  automated 
method  for  identification  of  equiaxed  a  in  SEM  micrographs. 

Following  these  two  initial  preprocessing  procedures,  the  stereological  procedures 
described  below  can  be  performed.  The  stereological  procedures  were  all  performed  on 
desktop  computers  running  Adobe  Photoshop  7.0  or  CS2  with  Fovea  Pro®,  a  commercial 
plug-in  set  available  from  Reindeer  Graphics,  Asheville,  NC.  Where  possible,  the 
following  procedures  were  fully  automated  using  the  action  and  droplet  functions 
available  in  Adobe  Photoshop. 

2.1  Stereological  Procedures  for  a+P  processed  Ti-AIIoys 

2.1.1  Mean  size  of  equiaxed  a 

In  addition  to  the  preprocessing  of  the  images  to  select  the  equiaxed  a,  it  was 
necessary  to  separate  the  individual  particles  from  clusters  of  the  equiaxed  phase  for  the 
calculation  of  its  mean  size  ( see  Figs.  4(c-d)).  This  additional  step  is  required  given  that  a 
cluster  of  equiaxed  a  particles  may  be  composed  of  regions  of  different  orientations,  as 
evidenced  by  the  variation  in  channeling  contrast,  or  have  internal,  high-angle  interfaces. 
Once  this  manual  step  is  completed,  the  mean  linear  intercept  was  determined  using  a 
Fovea  Pro™  plug-in.  This  plug-in  superimposes  a  grid  of  parallel  lines  on  a  binary 
image,  rotates  the  grid  by  successive  10°  increments  to  collect  measurements  in  multiple 
directions  (see  Fig.  5(b)),  and  records  both  the  mean  linear  intercept  and  the  mean  of  the 
inverse  of  the  intercepts  ( mean  inverse  intercept),  as  well  as  a  graphical  measure  of  the 
measurement  anisotropy  present  within  a  particular  micrograph.  It  is  important  to  make  a 
distinction  between  the  mean  linear  intercept  which  is  used  to  determine  the  mean  size  of 
the  equiaxed  a,  and  the  mean  inverse  intercept,  which  is  used  for  estimating  the  true 
thickness  of  Widmanstatten  a-laths.  The  mean  linear  intercept  provides  a  direct 
estimation,  with  a  shape  factor  where  appropriate,  for  geometrically  simple  shapes  with 
small  aspect  ratios  (e.g.,  spheres  and  cubes).  Conversely,  the  mean  inverse  intercept 
provides  for  an  estimation  of  the  true  three-dimensional  thickness  of  a  thin  membrane 
(e.g.,  oxide  layers  on  particles),  and  is  often  applied  to  the  estimation  of  thin,  layered 
structures  (e.g.,  pearlite,  thickness  of  Widmanstatten  a-laths).  Therefore,  the  mean  linear 
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intercept  is  appropriate  to  provide  an  estimation  of  the  equiaxed  a  present  in  these 
microstructures.  The  result  is  that  the  mean  size  of  the  equiaxed  a  particles  is 
approximately  equal  to  the  mean  linear  intercept  for  random  line  segments  in  the  a+p 
processed  Ti-6-4  microstructures.  Indeed,  for  true  spheres,  Dsphere  =  4/tx  dmean  [13]. 
However,  for  simplicity,  and  because  the  true  three-dimensional  shape  is  currently 
unknown,  the  mean  size  of  the  equiaxed  a  particles  is  not  multiplied  by  such  a  factor. 

2.1.2  Volume  Fraction  of  equiaxed  a  and  total  a 

After  the  preprocessing  of  the  high  resolution  images,  it  is  possible  to  use  those 
images  to  perform  an  area  fraction  measurement  to  determine  the  volume  fraction  of  the 
the  various  phases.  Although  there  are  two  separate  volume  fractions  of  interest  in  this 
microstructure,  that  of  the  equiaxed  a  phase  and  the  total  a  phase,  the  same  methodology 
can  be  used  to  determine  the  volume  fractions.  This  procedure  follows  the  classical 
stereological  procedures  where  the  point  fraction  (Pp),  line  fraction  (Ll),  area  fraction 
(Aa)  and  volume  fraction  are  equivalent.  Practically,  the  area  fraction  is  determined  using 
the  point  count  method,  where  the  high  resolution  SEM  micrographs,  such  as  the  one 
shown  in  Fig.  5(a),  are  essentially  high  density  point  grids,  where  every  pixel  is  a  point. 
This  method  can  be  automated,  as  was  done  for  the  analysis  of  the  volume  fraction  of  the 
equiaxed  a  (see  Fig.  5(c) )  and  the  total  volume  fraction  of  a  (see  Fig.  5(d)). 

2.1.3  Thickness  of  Widmanstatten  a-laths 

For  the  determination  of  the  thickness  of  the  Widmanstatten  a-laths  in 
transfonned  p,  it  is  essential  to  subtract  the  equiaxed  a  that  has  an  identical  contrast,  and 
therefore,  appears  the  same  in  the  binary  image.  A  Fovea  Pro™  plug-in  was  used  to 
perform  this  operation  (see  Fig.  5(e)).  Subsequently,  the  same  plug-in  used  to  measure 
the  size  of  the  equiaxed  a  particles  was  used  to  measure  the  mean  inverse  intercept 
length  of  the  a-laths.  As  described  in  section  2.1.1,  this  plug-in  superimposes  a  grid  of 
parallel  lines  on  a  binary  image,  rotates  the  grid  by  successive  10°  increments  to  measure 
in  all  directions  (see  Fig.  5(f)),  and  reports  both  the  mean  linear  intercept  and  the  mean  of 
the  inverse  of  the  intercepts.  Given  the  three-dimensional  shape  of  the  Widmanstatten  a- 
laths  (approaching  a  layered  structure),  the  mean  inverse  intercept  method  is  applied.  The 
relation  between  the  true  thickness  of  an  “infinite”  plate  of  finite  thickness  (e.g.,  a-laths) 
was  originally  derived  from  the  results  of  Gundersen  [14]  who  was  the  first  to  apply  the 
inverse  of  a  lineal  measurement  to  the  problem  of  a  thin  membrane  -  i.e.,  a  feature  finite 
in  one  dimension  and  infinite  in  two  dimensions.  In  this  work,  he  demonstrated  that 
while  the  probability  distribution  function  of  the  intercept  lengths  was  a  complicated 
asymptotic  shape,  the  probability  distribution  function  of  the  inverse  of  the  intercept 
lengths  was  a  triangle,  with  the  maximum  occurring  at  AT1  =  1  / ttnie.  Therefore,  as  the 
mean  (HX)mean  of  a  triangular  probability  distribution  function  is  2/3  of  the  max  (\/tlrue)  , 
the  maximum  can  be  estimated,  and  the  true  thickness  can  be  calculated  using  the 
following  relationship: 


thickness  = 


1 '5(V4 
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Work  now  suggests  that  the  probability  density  function  for  both  the  linear 
intercept  and  inverse  linear  intercept  for  the  case  of  infinite  plates  of  finite  thickness 
varies  from  the  original  work  of  Gundersen  [14]  ( see  Fig.  6  a-d).  Indeed,  as  can  be 
shown,  when  two  of  the  dimensions  become  finite  (i.e.,  the  case  of  a  one-dimension 
infinite  plate),  the  probability  density  functions  will  me  modified  to  such  a  degree  that  the 
mean  inverse  intercept  begins  to  deviate  from  the  value  of  0.67  (see  Fig.  7).  Although  the 
probability  distribution  function  of  the  inverse  intercept  no  longer  provides  a 
mathematically  simple,  tractable  solution,  the  detennination  of  the  mean  inverse  intercept 
((1  / X)mean)  is  still  currently  the  most  valid  approach  to  estimating  the  true  thickness  when 
the  ratio  of  the  two  dimensions  is  greater  than  -20: 1 .  It  is  expected  that  the  use  of  mean 
inverse  intercepts  will  break  down  for  any  case  where  the  feature  of  interest  does  not 
have  a  ratio  of  the  smallest  dimension  to  the  second  smallest  dimension  of  at  least  1 :20. 
While  these  observations  do  illustrate  limitations  with  the  mean  inverse  intercept  method, 
it  is  still  directly  applicable  for  the  case  of  Ti-based  alloys.  Additionally,  it  is  interesting 
to  note  that  the  calculated  mean  inverse  intercept  for  an  infinite  plate  is  0.637,  a  value 
that  is  less  than  5%  different  from  the  factor  of  0.667  derived  by  Gundersen.  It  is 
acknowledged  that  this  factor  will  change  with  particle  shape.  However,  for  consistency 
with  previous  work,  and  while  new  work  aimed  at  modeling  probability  distribution 
functions  for  particles  of  finite  shapes  is  currently  underway,  the  thickness  will  be 
estimated  by  the  above  relationship,  i.e.,  the  best  approximation  currently  available.  To 
the  best  knowledge  of  the  authors,  this  represents  the  first  time  where  Gundersen’ s  theory 
has  been  shown  to  be  not  applicable  to  all  thin-membrane  problems.  Additional  work  is 
currently  underway  to  develop  a  new  set  of  methods  for  thin  plates  of  various  sizes,  and 
will  be  the  subject  of  future  papers. 

3  Results  and  Discussion 

While  significant  variations  in  both  microstructure  and  properties  were  effected 
(Table  2),  the  primary  focus  of  this  paper  is  on  the  development  of  stereological 
procedures  that  can  be  used  to  describe  accurately  the  average  microstructural  features 
for  a+P  processed  a/p  Ti-alloys,  specifically  Ti-6A1-4V.  The  following  results  and 
discussion  will  focus  on  the  results  of  the  stereological  measurements,  with  particular 
emphasis  on  the  statistics  and  quality  of  the  measurements.  To  that  end,  four  figures  of 
merit  will  be  considered.  The  first  figure  of  merit  will  be  a  comparison  of  the  analysis  for 
each  of  the  two  identically  heat-treated  samples.  It  should  be  noted  that,  although 
microstructural  variation  between  the  two  samples  will  be  observed,  there  is  often  an 
associated  variation  in  the  properties,  as  little  as  ~1  MPa  ( 1  ksi)  and  as  great  as  42  MPa 
(6  ksi).  The  second  figure  of  merit  will  be  the  standard  deviation  for  the  four 
micrographs  that  were  recorded.  The  third  figure  of  merit  will  be  the  Coefficient  of  the 
standard  deviation  in  population  (VN/N,  %).  The  final  figure  of  merit  will  be  the 
coefficient  of  variation  (a/p  x  100,  %).  The  results  of  these  figures  of  merit  are  shown  in 
Table  3. 

Additionally,  select  microstructural  features,  as  described  by  their  mean,  will  be 
compared  with  both  processing  conditions  and  the  overall  yield  properties.  It  should  be 
noted  that  while  the  ultimate  goal  of  such  quantified  microstructural  features  is  the 
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population  of  a  database  to  be  used  for  the  training  and  testing  of  Bayesian  Neural 
Networks,  the  current  database  is  not  sufficiently  large  (i.e.,  only  17  sets  of  significant 
microstructural  variations)  to  develop  a  robust  neural  network. 

3. 1  Stereological  Methods 

3.1.1  Volume  Fraction  of  Phases 

In  addition  to  being  the  oldest  and  most  well  understood  stereological  procedure, 
the  detennination  of  the  volume  fraction  of  a  phase  can  also  prove  the  most  accurate 
under  certain  conditions.  The  primary  conditions  for  the  most  accurate  determination  of 
volume  fraction  are  that  the  plane  of  the  image  be  random,  and  that  the  presence  of  a 
particular  phase  is  equally  likely  at  all  points  within  a  microstructure.  Often,  phase 
transformations  satisfy  the  latter  condition.  Such  is  the  case  for  the  measurement  of  the 
total  volume  fraction  a  in  a+p  Ti-alloys.  However,  for  those  phases  that  preferentially 
form  at  certain  locations  of  the  microstructure,  the  measurement  of  volume  fraction  has 
as  many  uncertainities  as  other  standard  stereological  procedures.  For  example,  the 
presence  of  equiaxed  a  is  not  homogeneously  distributed  throughout  the  microstructure. 
Indeed,  the  conditions  required  for  the  formation  of  equiaxed  a,  namely  a  at  elevated 
temperatures  and  a  concurrent  strain,  are  not  homogeneously  distributed  throughout  the 
microstructure  (high  temperature  a)  or  sample  (s).  Therefore,  one  might  expect  low 
figures  of  merit  for  the  volume  fraction  equiaxed  a  as  measured  on  4  micrographs. 

The  quantified  microstructures  in  this  study  confirm  the  expected  outcome 
regarding  the  figures  of  merit  for  the  detennination  of  volume  fraction  of  both  total  a  and 
equiaxed  a.  The  figures  of  merit  indicating  the  data  of  the  higest  fidelity,  including 
notably  the  lowest  sample-to-sample  variation  (1.17%)  and  lowest  coefficient  of  variation 
(0.98%),  were  those  for  the  measurements  of  the  volume  fraction  total  a.  Conversely, 
the  lowest  figures  of  merit,  including  both  the  greatest  sample-to-sample  variation 
(1 1.2%)  and  greatest  coefficient  of  variation  (9.56%)  were  those  for  the  measurements  of 
the  volume  fraction  equiaxed  a.  As  described  above,  the  heterogeneous  nature  of 
equiaxed  a  particle  fonnation  is  the  cause  for  the  uncertainity.  Therefore,  it  is  expected 
that  the  figures  of  merit  would  improve  with  an  increase  in  the  sampling  statistics  (e.g., 
more  micrographs)  and  the  corresponding  decrease  in  sample  standard  deviations. 

3.1.2  Mean  Equiaxed  a  particle  size 

The  stereological  procedure  for  the  determination  of  the  mean  equiaxed  a  particle 
size  results  in  reasonable  figures  of  merit,  as  shown  in  Table  3.  It  is  interesting  to  note 
that  the  three  principal  figures  of  merit  (sample-to-sample  variation  coefficient,  standard 
deviation  of  the  population,  and  average  coefficient  of  variation)  are  similar  in  magnitude 
(5.8%,  5.7%,  and  4.8%,  respectively).  These  are  the  second  best  figures  of  merit,  after 
the  volume  fraction  of  total  a.  It  is  expected  that  each  of  the  figures  of  merit  would 
decrease  with  an  increase  in  sample  population  from  the  current  population  of  320-600 
particles  counted  per  sample  (80-150  particles/micrograph,  4  micrographs).  However, 
10-20  micrographs  would  be  required  per  sample  to  increase  the  total  population  to 
approximately  1600  particles,  resulting  in  a  decrease  in  the  standard  deviation  of  the 
population  to  +  2.5%.  Such  sampling  may  or  may  not  be  required. 
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3.1.3  Thickness  of  Widmanstatten  a-laths 

Based  on  the  figures  of  merit  in  Table  3,  the  method  to  detennine  the  thickness  of 
the  Widmanstatten  a-laths  is  a  relatively  robust  method.  Indeed,  it  is  suggested  that  the 
results  obtained  from  4  micrographs  of  an  a+p  processed  Ti  microstructure  are  more 
accurate  than  those  from  a  p-processed  Ti  microstructure.  This  expectation  is  the  result  of 
an  improvement  in  the  random  sampling  of  different  orientations  that  a+p  processed  Ti- 
alloys  afford,  given  the  smaller  regions  of  P  regions  of  different  orientations  that  exist 
prior  to  transformation.  The  figures  of  merit  are  misleading,  given  that  the  resolution  of 
these  micrographs  is  0.05  pm  /  pixel.  Indeed,  an  adjustment  of  a  single  pixel  to  an 
average  lath  thickness  (2/3*pixel)  would  result  in  a  value  that  is  identical  to  the  sample- 
to-sample  variation  and  average  standard  deviation  that  have  been  experimentally 
observed,  0.03  pm  and  0.024  pm,  respectively.  By  extension,  for  the  thickness  of  the 
Widmanstatten  a-laths,  it  is  expected  that  the  coefficients  of  variations  can  not  be  less 
than  ~5-7%,  as  these  are  limited  by  the  resolution  of  the  microscope.  Therefore,  while  it 
is  acknowledged  that  the  actual  stereological  pre-factor  might  change  the  accuracy  of  the 
calculations,  the  precision  is  only  limited  by  the  resolution  in  the  micrographs.  Similarly, 
it  is  important  to  compare  the  precision  afforded  by  the  high  resolution  SEM  micrographs 
(500x  in  this  case)  to  those  of  optical  micrographs.  Although  the  highest  pixel  resolution 
in  optical  microscopes  is  often  0.1  pm,  the  highest  practical  resolution  (0.2  pm)  is 
governed  by  the  Rayleigh  Criterion  (the  criterion  for  the  minimum  resolvable  detail  in  an 
image)  and  occurs  for  an  optimized  optical  microscope  with  a  light  wavelength  of  550 
mn  and  a  numerical  aperture  of  1.35  to  1.40.  Clearly,  the  uncertainty  afforded  by  using 
electron  microscopy  is  significantly  lower  than  what  would  be  afforded  using  optical 
microscopy. 

3.2  Correlation  with  processing 

One  benefit  to  the  development  of  such  databases  is  that  it  affords  the  possibility 
of  investigating  a  variety  of  interdependencies,  including,  but  not  limited  to,  the  influence 
of  microstructure  on  properties  and  the  influence  of  processing  on  microstructure.  While 
the  variation  in  thermomechanical  processing  is  significant  in  this  database,  it  is  possible 
to  keep  most  processing  variables  constant  and  investigate  the  influence  of  one 
processing  variable  on  microstructure.  For  example,  while  one  would  expect  the  solution 
temperature  and  cooling  rates  to  influence  the  size  of  the  equiaxed  a  particles  and 
Widmanstatten  a-laths,  it  is  beneficial  to  exercise  the  databases  to  demonstrate  this.  Figs. 
8  (a-c)  show  the  influence  of  isolated  processing  variables  on  such  microstructural 
features.  Fig.  8(a)  shows  that  as  the  cooling  rate  from  an  aging  temperature  increases 
(both  593°C  and  704°C  shown),  the  thickness  of  the  Widmanstatten  a-laths  decreases. 
Fig.  8(b)  shows  that  as  the  solution  temperature  decreases,  the  volume  fraction  of 
equiaxed  a  increases  and  Fig.  8(c)  shows  that  as  the  cooling  rate  from  two  separate 
temperatures  decreases,  the  volume  fraction  of  equiaxed  a  increases.  It  is  noted  that 
these  plots  show  the  average  and  standard  deviations  of  the  stereological  measurements 
for  each  of  two  identically  processed  samples.  Thus,  two  discrete  average  data  points  and 
their  uncertainties  are  provided  for  each  value  of  a  particular  processing  variable. 
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Interestingly,  the  relative  uncertainties  of  the  Widmanstatten  a-laths  and  equiaxed  a- 
particles  are  clearly  observed  when  comparing  Figs.  7(a)  and  7(c),  with  the 
Widmanstatten  a-laths  exhibiting  a  lower  sample-to-sample  variability  and  a  lower 
uncertainty.  These  two  trends  are  expected  given  the  current  understanding  and  models 
for  the  nucleation  and  growth  of  the  equiaxed  a  particles. 

3.3  Correlation  with  properties 

Figs.  9  (a-d)  show  the  raw  data  correlations  between  the  ultimate  tensile  strength 
and  a  particular  microstructural  variable.  As  can  be  observed,  while  it  might  be  tempting 
to  conclude  that  an  increase  in  either  equiaxed  a  particle  size  or  Widmanstatten  a-lath 
width  results  in  a  decrease  in  strength,  the  scatter  in  the  data  indicates  that  such  trend 
plots  do  not  allow  a  clear,  unambiguous  interpretation  of  the  influence  of  a  single 
microstructural  variable  and  the  resultant  property,  in  this  case,  ultimate  strength.  The 
scatter  in  the  plots  is  not  an  indication  of  the  quality  of  the  quantified  microstructural 
feature.  Rather,  it  is  a  direct  consequence  of  the  microstructural  interdependencies.  For 
example,  consider  the  five  samples  highlighted  in  Table  4.  Each  sample  has  an 
equivalent  ultimate  tensile  strength  (within  7  MPa  or  1  ksi),  yet  each  of  the  four 
microstructural  variables  changes  in  a  complicated,  interrelated  fashion.  Such 
interdependencies  require  more  robust  modeling  schemes,  such  as  Bayesian  Neural 
Networks,  which  require  a  larger  database  than  that  presented  in  this  paper.  These  are  the 
subject  of  ongoing  work. 

4  Conclusions 

Standardized  stereology  procedures  have  been  developed  for  the  quantification  of 
microstructural  features  in  a+[!  processed  Ti-alloys,  including  the  total  volume  fraction 
alpha,  the  volume  fraction  of  equiaxed  alpha,  the  mean  size  of  the  equiaxed  a  particles, 
and  an  estimation  for  the  thickness  of  the  Widmanstatten  a-laths  in  the  transformed  [! 
region.  The  development  of  these  procedures,  and  the  subsequent  automation  or  semi¬ 
automation  through  software  such  as  Adobe  Photoshop™  and  Fovea  Pro™,  have 
decreased  or  eliminated  the  bias  (and  associated  errors)  that  occurs  in  more  classical, 
fully  manual  stereology  procedures.  The  uncertainities  associated  with  each  of  the  four 
procedures  have  been  discussed. 
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Figure  Captions: 

Fig.  1:  (a)  basketweave  microstructures  from  p-processed  Ti-6A1-4V  and  (b)  equiaxed 
microstructures  from  a+p  processed  Ti-6A1-4V. 
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Fig.  2:  (a)  Example  micrograph  (b)  histogram  from  alpha  region  with  channeling 

contrast  (c)  histogram  from  alpha  region  without  channeling  contrast  (d)  histogram  from 
complete  microstructure  and  (e)  threshheld  image 

Fig.  3:  (a)  small  equiaxed  a,  low  volume  fraction  (b)  large  equiaxed  a,  small  volume 
fraction  (c)  large  equiaxed  a,  large  volume  fraction 

Fig.  4:  (a)  Original  backscattered  electron  micrograph  (b)  original  automated 
identification  of  equiaxed  a  illustrating  over-identification  (c)  manual  refinement  of 
automated  identification  and  (d)  manual  separation  of  clusters  of  equiaxed  a  particles 

Fig.  5:  Original  image  (a)  and  representative  processed  images  for  the  measurement  of 
(b)  equiaxed  a  particle  size  (c)  volume  fraction  of  equiaxed  a  (d)  total  volume  fraction  of 
a  and  (e)  isolation  and  (f)  thickness  measurement  of  the  Widmanstatten  a-laths. 

Fig.  6:  Probability  distribution  functions  vs  the  intercept  length  for  (a)  Gundersen’s 
derivation  of  membranes  (b)  distribution  through  an  infinite  plate  of  finite  thickness  (c) 
comparison  of  the  two;  and  probability  distribution  functions  vs.  the  inverse  intercept 
length  for  (d)  Gundersen’s  derivation  of  membranes  (e)  infinite  plates  of  finite  thickness 
and  (f)  a  comparison  of  both. 

Fig.  7:  Calculated  mean  inverse  intercept  from  the  probability  density  function  of  inverse 
intercept  for  1-D  infinite  plates  where  the  two  other  dimensions  (width  and  thickness)  are 
given  by  their  ratio. 

Fig.  8:  Influence  of  processing  variables  on  selected  microstructural  features,  including 
(a)  Widmanstatten  a-lath  thickness  vs.  cooling  rate  from  the  aging  temperature  for 
samples  subjected  to  a  total  strain  of  80%  (b)  solution  temperature  vs.  the  volume 
fraction  of  equiaxed  a  for  samples  subjected  to  a  total  strain  of  40%  and  cooled  at 
3.7°C/sec  and  (c)  the  volume  fraction  of  equiaxed  a  vs.  the  cooling  rate  from  the  aging 
temperature  for  samples  subjected  to  a  total  strain  of  80%. 

Fig.  9:  Raw  data  scatter  plots  for  the  UTS  vs.  the  four  microstructural  variables  quantified 
in  a+p  processed  a/p  Ti  alloys  (here  Ti-6A1-4V). 

List  of  Tables: 

Table  1:  Full  details  of  thermomechanical  processing. 

Table  2:  Ranges  of  microstructural  variables  as  determined  using  stereology  and  the 
ranges  of  properties  in  the  developed  database. 

Table  3:  Figures  of  Merit. 

Table  4:  Samples  with  similar  properties  but  different  microstructural  features. 


11 


Fig.  1:  (a)  basketweave  microstructures  from  p-processed  Ti-6A1-4V  and  (b)  equiaxed 
microstructures  from  a+p  processed  Ti-6A1-4V. 
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Fig.  2:  (a)  Example  micrograph  (b)  histogram  from  alpha  region  with  channeling 

contrast  (c)  histogram  from  alpha  region  without  channeling  contrast  (d)  histogram  from 
complete  microstructure  and  (e)  threshheld  image 
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Fig.  3:  (a)  small  equiaxed  a,  low  volume  fraction  (b)  large  equiaxed  a,  small  volume 
fraction  (c)  large  equiaxed  a,  large  volume  fraction 
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Fig.  4:  (a)  Original  backscattered  electron  micrograph  (b)  original  automated 
identification  of  equiaxed  a  illustrating  over-identification  (c)  manual  refinement  of 
automated  identification  and  (d)  manual  separation  of  clusters  of  equiaxed  a  particles 
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Fig.  5:  Original  image  (a)  and  representative  processed  images  for  the  measurement  of 
(b)  equiaxed  a  particle  size  (c)  volume  fraction  of  equiaxed  a  (d)  total  volume  fraction  of 
a  and  (e)  isolation  and  (f)  thickness  measurement  of  the  Widmanstatten  a-laths. 
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Gundersen's  Approximation 
for  membranes  covering  particles 


(b) 


Distribution  of  intercept  lengths  for 
random  lines  through  a  membrane 
with  infinitely  parallel  interfaces 


(c) 


Distribution  of  intercept  lengths  for 
random  lines  through  a  membrane 
with  infinitely  parallel  interfaces 


(a) 


(d) 


Intercept  length  ( >.) 

Distribution  of  inverse  intercept  lengths 
as  derived  by  Gundersen 


(e) 


Intercept  length  ( X) 


Distribution  of  inverse  intercept  lengths 
of  random  lines  through  an  infinite  plate 
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Fig.  6:  Probability  distribution  functions  vs  the  intercept  length  for  (a)  Gundersen’s 
derivation  of  membranes  (b)  distribution  through  an  infinite  plate  of  finite  thickness  (c) 
comparison  of  the  two;  and  probability  distribution  functions  vs.  the  inverse  intercept 
length  for  (d)  Gundersen’s  derivation  of  membranes  (e)  infinite  plates  of  finite  thickness 
and  (f)  a  comparison  of  both. 
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Mean  Inverse  Intercepts  for  1-D  infinite  plate 
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Fig.  7:  Calculated  mean  inverse  intercept  from  the  probability  density  function  of  inverse 
intercept  for  1-D  infinite  plates  where  the  two  other  dimensions  (width  and  thickness)  are 
given  by  their  ratio. 
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(a) 


Influence  of  cooling  rale  on 
thickness  of  Widmanstatten u-laths 


Influence  of  Solution  Temperature 
on  volume  fraction  equiaxed  a 


Influence  of  cooling  rate  on 
volume  fraction  equiaxed  << 
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Fig.  8:  Influence  of  processing  variables  on  selected  microstructural  features,  including 
(a)  Widmanstatten  a-lath  thickness  vs.  cooling  rate  from  the  aging  temperature  for 
samples  subjected  to  a  total  strain  of  80%  (b)  solution  temperature  vs.  the  volume 
fraction  of  equiaxed  a  for  samples  subjected  to  a  total  strain  of  40%  and  cooled  at 
3.7°C/sec  and  (c)  the  volume  fraction  of  equiaxed  a  vs.  the  cooling  rate  from  the  aging 
temperature  for  samples  subjected  to  a  total  strain  of  80%. 
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Volume  Fraction  Equiaxed  a  Volume  Fraction  Total  alpha 


Fig.  9:  Raw  data  scatter  plots  for  the  UTS  vs.  the  four  microstructural  variables  quantified 
in  a+p  processed  a/p  Ti  alloys  (here  Ti-6A1-4V). 
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